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The first-order superconducting fluctuation corrections to the thermal conductivity of a granular 
metal are calculated. A suppression of thermal conductivity proportional to T c j (T — T c ) is observed 
in a region not too close to the critical temperature T c . As T ~ T c , a saturation of the correc- 
tion is found, and its sign depends on the ratio between the barrier transparency and the critical 
temperature. In both regimes, the Wiedemann-Franz law is violated. 



I. INTRODUCTION 



In normal metals, in the presence of BCS interaction, 
electrons can form Cooper pairs even for temperatures 
T larger than the critical temperature T c . As T > T c , 
the pairs have a finite lifetime, the Ginzburg-Landau 
(GL) time, inversely proportional to the distance from 
the critical temperature tql ~ (T-T c ) . These super- 
conducting fluctuations strongly affect both the thermo- 
dynamic and transport properties and since many years 
they are widely studied both theoretically and experi- 
mentally. 

The first analysis of fluctuation corrections has been 
performed on electrical conductivity where the pair- 
ing leads to three distinct contributions named the 
Aslamazov-Larkin (AL), the Maki-Thompson (MT) and 
Density of States (DOS) terms. In the first one, the for- 
mation of Cooper pair leads to a parallel superconduct- 
ing channel in the normal phase; the second takes into 
account the coherent scattering off impurities of the (in- 
teracting) electrons; finally, the third one is due to the 
rearrangement of the states close to the Fermi energy 
since electrons involved in pair transport are no longer 
available for single particle transport. Both the AL and 
MT terms lead to an enhancement of the conductivity 
above T c , on the contrary, the DOS correction is of op- 
posite sign. 

The analysis of superconducting fluctuation correc- 
tions to thermal conductivity dates back to the early 
1960s, when Schmid 2 and Caroli and Maki 3 found an 



expression for the heat current in the framework of the 
phenomenological time dependent GL theory, (TDGL). 
More recently, a complete analysis was performed, in the 
same framework of the TDGL, by Ussishkini. Abrahams 
et aim' first pointed out the divergence of the thermal con- 
ductivity in the vicinity of the critical temperature due to 
the opening of the fluctuation pseudogap in the density 
of states (DOS) energy dependence in the homogeneous 
case. Niven and Smith have shown& that Abrahams's 
DOS correction [« G£ln(l/e), e = (T-T c )/T c , Gi being 
the so-called Ginzburg-Levanyuk parameter] is exactly 
compensated by the regular Maki-Thompson (MT) one; 
hence, all singular first order fluctuation corrections are 
cancelled out. The only surviving contribution to heat 
conductivity, the Aslamazov-Larkin (AL) one, is non- 
singular in temperature. Therefore, in bulk metals, no 
singular behaviour of the heat current is expected at the 
metal-superconductor phase transition. 

In this paper we are interested in the superconduct- 
ing fluctuation corrections to the thermal conductivity 
in a granular superconductors, an ensemble of metallic 
grains embedded in an insulating amorphous matrix and 
undergoing a metal-superconductor phase transition due 
to the existence of pairing interaction inside each grain. 
The electrons can diffuse in the system due to tunneling 
between the grains. Experimentally, this kind of systems 
have been investigated, for example, in Ref. 0. Each Al 
grain has an average dimension of 120A, while the sample 
has a linear dimension of the order of mm, that is, much 
larger than the superconducting coherence length. The 
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reason for studying thermal transport in granular metals 
is that, depending on the temperature regime, a radically 
different behaviour, as compared with the homogeneous 
case, may emerge. In fact, in granular material (a simi- 
lar situation occurs in layered superconductors) the AL 
and MT contributions are of higher order in the tunnel- 
ing amplitude as compared to the DOS. This effect has 
been observed, for example, in the electrical^ and the 
optical conductivity 9 ! of layered superconductors and the 
electrical conductivity^ of granular systems. Indeed, in 
granular superconductors there is a temperature region in 
which a singular correction due to superconducting fluc- 
tuations for a quasi-zero-dimensional system dominates 
the behaviour of the thermal conductivity; such a correc- 
tion can be either negative or positive, depending on the 
ratio between the barrier transparency and the critical 
temperature T c . When the temperature approaches T c , 
the behaviour observed in homogeneous systems is recov- 
ered, and the divergence will be cut off to crossover to 
the regular behaviour. Moreover, a significant difference 
with the homogeneous systems is present, the constant 
correction at T = T c being either negative or positive de- 
pending on the above-mentioned ratio. For some choice 
of the parameter, a non-monotonic temperature depen- 
dent behaviour of the correction is possible. 

A phenomenological approach to granular supercon- 
ductors has been proposed long agoii*i£, while the micro- 
scopic theory has only been formulated very recentlyiiS. 
The difference between bulk and granular microscopic 
theory is mainly based on the renormalization of the su- 
perconducting fluctuation propagator due to the pres- 
ence of tunneling. This renormalization accounts for 
the possibility that each electron forming the fluctuat- 
ing Cooper pair tunnels between neighbor grains during 
the Ginzburg-Landau time. 

The paper is organized as follows. In Section[nl we de- 
scribe and formulate the model. Section IlIII contains the 
main steps and assumptions of the calculation of fluctua- 
tion propagator. Its expression, calculated in Ref. lid is 
given explicitly at every order in tunneling in the ladder 
approximation. By means of that, DOS, MT, and AL 
corrections are evaluated. For each of those corrections, 
an explicit form for the response function is presented. In 
the final section, we discuss the overall behaviour of the 
fluctuation corrections to thermal conductivity as a func- 
tion of temperature. For temperatures sufficiently far 
from T c , the system behaves as in the zero-dimensional 
case. In this region, the correction to the heat conductiv- 
ity has a singular behaviour: \5k\ oc Kq/ (grt) , where kq 
is the classical Drude conductivity for a granular metal, 
and it reads 

= y 9 T a 2 - d T , (1) 

a being the size of a single grain, d the dimensionality of 
the system and e = (T — T c ) /T c the reduced tempera- 
ture. We defined the dimensionless macroscopic tunnel- 
ing conductance = [{it /2)tvp\ 2 , with Vp the electronic 
density of states at the Fermi level, and t the hopping 



energy. On the other hand, when the correlation length 
increases until the distance between two nearest neighbor 
grains, the tunneling becomes important and the correc- 
tion, exactly at the critical temperature, reduces to a 
constant 

5«=JLf±^-8V. (2) 
zgr \2tt T c tt z J 

Connections with the homogeneous metal results are dis- 
cussed. In the appendix, we briefly review the evaluation 
of the superconducting fluctuation propagator in a gran- 
ular metal. Throughout the paper, we set h — kp = 1. 

II. THE MODEL 

We consider a d— dimensional array of metallic grains 
embedded in an insulating amorphous matrix, with im- 
purities on the surface and inside each grain. Even if 
the analytical model we use is for a perfectly ordered d- 
dimensional matrix, the results we found still hold for an 
amorphous one. Indeed, one can imagine different pos- 
sible configurations of spatial position of grains in the 
lattice, that is, different disordered configurations. Con- 
sequentely, the hopping matrix shall vary for each sam- 
ple. By performing the average over disorder, one gets a 
model with the same value of coordination number and 
hopping energy, t, for different configurations. In other 
words, our description is correct until the system can be 
described by a dimensionless tunneling conductance, gx, 
on a scale which is much bigger than the typical linear 
dimension of the grains, a, but smaller than the macro- 
scopic dimension of the whole sample. 

The Hamiltonian of the system reads 

H = H + H P + H T • (3) 

Hq and Hp describe the free electron gas and the pairing 
Hamiltonian inside each grain, respectively 

Hq = £ i,ka i k Cli t k + -Himp , (4) 

Hp = -X ^2 a\ k a\ _ k a,i -ya^ , (5) 

i,kk f 

where i is the grain index, and d\ k (a.j,fc) stands for cre- 
ation (annihilation) operator of an electron in the state 
k = (k, f) or — k = (— k, J.). The term Hi mp describes the 
electron elastic scattering with impurities. The interac- 
tion term in Eq. Q contains only diagonal termsi^. Such 
a description is correct in the limit 

6 < A < E T , (6) 

where S ~ Vp 1 is the mean level spacing and the smallest 
energy scale in the problem, and A the (BCS) supercon- 
ducting gap of a single grain, supposed equal for each 
of them. Et — D/a 2 is the Thouless energy, D being 
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FIG. 1: Diagram for the thermal conductivity in granular 
metals. The solid lines are impurity-averaged single-electron 
Green's functions with the specified momentum and Matsub- 
ara's frequency, and belonging to the grain i and j. The 
vertices are discussed in the text. 



are shorthand notations for e„ + (i/2r)signe„ and e n +w„, 
respectively. In the latter equation, we considered the 
tunneling amplitude uniform and momentum indepen- 
dent, t p J = t. The thermal conductivity for free electrons, 
KOi Eq. |T|I. is given by the diagram in Fig. ^ where, 
as usual, Green's function is G(e„,p) = l/[is n — £(p)], 
and each vertex contributes as i2at(e n + u> u /2). Elec- 
trical conductivity reads Co = e 2 (8/ir)gTa 2 ~ d ; therefore, 
the Lorenz number is Lq = kq/<jqT = 7r 2 /3e 2 . 

III. SUPERCONDUCTING FLUCTUATION 
CORRECTIONS TO THERMAL 
CONDUCTIVITY 



the intragrain diffusion constant. Under the previous as- 
sumption, Eq.©, one can safely neglect off-diagonal 1/g 
corrections, where g is the dimensionless conductance of 
a grain, g = Et/S. Equation © is equivalent to the 
condition a <C £o, where £o = y/D/T c is the dirty super- 
conducting coherence length; then, Eq.© describes an 
ensemble of zero-dimensional grains. In addition, Eq. © 
states that the energy scale, r _1 , with r being the mean 
free time, related to -ff; mp is much larger than A. 

The grains are coupled by tunneling. The tunneling 
Hamiltonian is written as (t <C Et) 



EE [ 



t pq a) a 



H.C 



(7) 



We assume that the momentum of an electron is com- 
pletely randomized after the tunneling. Finally, assuming 
that the system is macroscopically a good metal, t >■ 5, 
we can safely neglect the Coulomb interaction, it being 
well screened^, and weak localization corrections too, at 
least for not too low temperatures 17 , i.e. when T < grS. 

The tunneling heat current operator is given as 



j(hcat) = ia 



EE[ 

pq<T 



e t pq a} a 



H.c. 



(8) 



where e n is the Matsubara frequency of the electron in- 
volved in the transport. 

In linear response theory, the heat conductivity is de- 
fined as 



lim 



to„T 



(9) 



-i0+ 



where Q( hcat ) (wj) is the linear response function to an 
applied temperature gradient: 



)(hcat) 



n-\-i> 



X j(dp)G (i n+u , p) J(dq)G(£ n ,q) , (10) 

where G(e n ,p) is the exact Matsubara Green's function 
of an electron in a grain, (dp) — [d d p/ (2ir) d ], i n and £ n + v 



At temperatures above but not far from the critical 
one, superconducting fluctuations allow the creation of 
Cooper pairs that strongly affect transport. In other 
words, fluctuations open a new transport channel, the 
so-called Cooper pair fluctuation propagator, Ref. Q. It is 
such a contribution that gives rise to corrections to both 
the electrical and thermal conductivity. 

With respect to the bulk case, the propagator is renor- 
malized by the tunneling, and as explained in the ap- 
pendix, it takes into account the possibility that each 
electron forming the cooper pair can tunnel from one 
grain to another, without loosing the coherence. 

The expression for the superconducting fluctuation 
propagator for a granular metal, calculated in Ref. [ToL 



A K = - 



1 



TrjfVj 
8T C 



^(1-7k) : 



(11) 



where K is the wave vector associated with the lattice 
of the grains, fi M is a bosonic Matsubara's frequency, 
and z the number of nearest neighbor grains. The func- 
tion 7k = (l/^) 5Z a e lK a is the so-called lattice structure 
factor, where a is a vector connecting nearest neighbor 
grains. The main steps of the calculation of Eq. Ullfl . 
done in Ref. ITol are reviewed for completeness in the 
appendix. 

The various contributions to thermal conductivity are 
shown in Figs. [2] and [31 

The correction due to the density of states renormal- 
ization, Fig. Ufa), is the only one which is present 
even in absence of tunneling; therefore, for temperatures 
T — T c ^> gr5, we expect this term to give a significant 
contribution to the thermal conductivity. For lower tem- 
peratures, the bulk behaviour will be recovered. 

The MT correction, represented in Fig. |2fb), can be 
evaluated using the same procedure as in the case of the 
DOS one. It is important to stress that the sign of linear 
response function is the same as for the DOS: in fact, 
the energies of electrons entering the diagram from op- 
posite sides have opposite signs but the same happens to 
their velocities. In the case of electrical conductivity, the 
sign of linear response function is opposite. It is this dif- 
ference that ultimately results in the cancellation of two 



4 




and 



FIG. 2: Diagrams representing various fluctuation contribu- 
tions to the thermal conductivity, (a) density of states con- 
tribution, (b) Maki-Thompson contribution. The solid lines 
are impurity-averaged single-electron Green's functions, wavy 
lines represent the fluctuation propagator and the shadowed 
areas are Cooperon vertex corrections. Crossed circles repre- 
sent tunneling vertices. 



identical contributions in the thermal conductivity 6 . 

Let us finally comment on the AL contribution, given 
by the diagrams in Fig. [3] It is well known, in the case of 
homogeneous metals, that such a correction to the ther- 
mal conductivity is not singulas 6 ^. We will show briefly 
that in the case of granular metals this correction van- 
ishes^ in the static limit too, but not in the dynamical 
one, giving an important and characteristic contribution 
to the total correction. 

In the following paragraphs, we present the evaluation 
of corrections to thermal conductivity due to different 
diagrams. 



A. Density of states correction 

The diagram for the DOS correction is given in Fig. 
Ufa) and the corresponding response function can be writ- 
ten as 

Q(DOS) {uJv) = T 2 t 2 a 2 J- Ay £ (%, Uu) , (12) 

(hj) ^ 

where 

11(0,^, U! v ) = ^""^ A 2 (gK+jy, S-n-v+ii) i^n + £n+v) 

xI(e n ,Q IJ ,,( J j u ) , (13) 



I (e n ,fl,i,u u ) = J (dp) Gq (p, e n+v ) G (p, e_ n _„ +At ) 



x / (dp')G (p',e n ). 



(14) 



We introduced the Cooperon vertex correction, 
A (£1,22) = l/ r l £ i~ £ 2| m the zero-dimensional 
limit and without tunneling corrections 2 ^. The main 
contribution to singular behaviour comes from "classi- 
cal" frequencies, |f2 M | <C T c : consequently, we will take 
the so-called static limit, fl^ = 0, in the calculation of 
correction. This will be true also for the Maki-Thompson 
correction in the next paragraph. In the dirty limit, we 
can neglect all the energy scales in the electronic Green's 
function in comparison with 1/r S> T, and the factor 
/ (s n , 0, ui v ) turns out to be 

I (e„, 0, w„) = -2 (■kv f t) 2 [9 (e n e n+u ) - 9 (-e n e n+v )] . 

Inserting the previous expression in Eq. (|13|) , we are left 
with the sum over the electronic Matsubara frequencies. 
It is straightforward to check that the only contribution 
linear in w v is given by E (0, w u ) — —w v %v%. By means of 
Ea.Q12|). we obtain the general form for the DOS response 
function after the analytical continuation 



q(DoS) / 



-ioj) -g T Ta A 2 

7T ^ — ' 



(0) , (15) 



where we also took into account the multiplicity of the 
DOS diagrams. The corresponding correction to heat 
conductivity is given by 



Sk (DoS) 



K 



[ {dK) i 



.(16) 




We took the lattice Fourier transform and defined 
the reduced temperature e = ln(T/T c ) ~ (T - T c )/T c . 
(dK) = [a d / (2n) d ]d d K is the dimensionless measure of 
the first Brillouin zone. Close to T c , the integral takes its 
main contribution from the small momentum region and 
we recover the bulk DOS behaviour as 



(17) 



B. Maki-Thompson correction 

The MT correction, [Fig. reads 

Q< MT > (uj v ) = a 2 Tt 2 E A v ^ E "») > (18) 

where 

) A (e n , £- n +/i) 
YI{e n ,%,u3 v ) , (19) 
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FIG. 3: Diagrams of the blocks appearing in the Aslamazov- 
Larkin contribution to thermal conductivity, (b) has a double 
molteplicity, since the bare tunnel vertex can stay on both side 
of the block. 



and 

I(e n ,fl^,LJ 1/ ) = J (dp) G (p, £„+„) G (p, E-n-u+fj) 

x / (dp') G (p',e n ) G (p', £_„+„) (20) 



Using the same procedure outlined above to calculate the 
DOS correction, we get 



Sk^ 3 1 g T S 



Kq 



tt 2 9T T c J BZ 



(dK) 



7K 



e + z^(l-7 K ) 



(21) 



As expected, the MT correction has the same singular 
behaviour as the DOS but opposite sign. On the other 
hand, because such a correction involves the coherent 
tunneling of the fluctuating Cooper pair from one site 
to the nearest neighbor, it is proportional to the lattice 
structure factor 7k: due to this proportionality, in the 
regime T ~ T c 3> gT$, the correction vanishes because 
J BZ (dK)jK = 0. Let us stress again that this is not the 
case for the DOS correction, which in this regime behaves 
as -(l/g)(E T /T c )(l/e). 



C. Aslamazov-Larkin correction 

The AL diagrams can be built up by means of blocks 
in Fig. |3 by considering all their possible combinations 
in pairs. For a sake of simplicity, we will call the first 
block, Fig. Ota), -Bi, and the second one B-i- Finally, one 
has three different kind of diagrams: the first one, with 
two -Bi-type blocks; the second one with two i?2-type 
blocks, and the latter, with both of them. Because of the 
double molteplicity of i?2-type block, one has a total of 
nine diagrams contributing to thermal conductivity. In 
the following, first we evaluate the analytical expression 
of B\ and Bi in the static approximation, then in the 
dynamical one, giving the expression of the total AL cor- 
rection. 

The general expression of response function for the AL 



diagrams reads 

q( al > K) = TVt 4 ]T A y (*W») A ™< ( n *) 

(M) f2 >. 

x.Bi e f t (uj v , Qfj) -B r i g ht (u) v , O^) , (22) 

where Bi f t and -B r i g ht can be either B\ or i?2-type. 
-Bi block reads 

n £n+fc') A (^"7t+i/i £/4 — n) A (&ni £/i— n) 

x J(dp)G (p,e n+J/ )G (p,e M _ n ) 



x / (dp')Go (p', £»-„) G (p',s„) . 



(23) 



Taking the integrals over the Fermi surface, in the static 
approximation, we get 

Bi(cj v ,0) = (2TTv F T) 2 ^26(s n+u s n )(e n + e n+u ) 
x A (e n +„, — e n ) A (e„, — e n ) 



(2nv F ) 2 

E +E 

_e n <—u) v e n >0_ 



~-n-\-i> 



\e n +u + £n\ |2e„| 
manipulating the sum, it is easy to see that 

BiK,0) = (2irv F ) 2 V -L 



(24) 



0<e n <u, 



= (2-KVpf 



r \ 2ttT 2 



(25) 



In the same way as sketched above, one can show, al- 
ways in the static approximation, that the block B2 van- 
ishes identically. Then, all the diagrams containing B 2 - 
type blocks do not give any contribution. Since the only 
AL diagram with two Bi-type block is proportional to 
the square of Eq. I|25|) , it is quadratic in the external fre- 
quency oj, and therefore vanishes identically in the limit 
lu -> 0. 

To evaluate the first non vanishing AL correction, one 
has to consider the dynamical contribution. In such a 
case, the B2 block, for instance, reads 

S 2 (o) v ,fl M ) = (e n + £„+„) A (e 

x J ' (dp')Go (p', e^-n) G (p', £„) G (p', £„+„) 

(26) 



x / (dp)G (p,e n+v ) • 

In the evaluation of the block, because of the pole 
structure of fluctuation propagator, one can neglect the 
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u) v dependence^, and keep just the one in fi^. The cal- 
culation of the integrals and the sums in the latter equa- 
tion is, in the dynamical approximation, a little bit more 
cumbersome. One has to take into account the different 
possible signs of and e n . Finally, Eq. it2T)|l reads 



B 2 (0,^) = -2(tt^) 2 ^ 



2e, 



(2e„ - 



+ %«))] (27) 



9{x) being the step function. 

By taking the lowest order in the bosonic frequency 
flu, one gets the result for the block 



B 2 (0, fi. 



2 V 2T J ' 



(28) 



In the same way, one can evaluate also B\ with the 
result 



B\ (0,O M ) = -2 J B 2 (o,n M ) 



(29) 



which is consistent with the homogeneous caseA^. The 
sum over fi^ in the response function can be performed 
by writing the sum as an integral^, and exploiting the 
properties of the pair correlators. 

Finally, the AL dynamical correction to thermal con- 
ductivity reads 



Sk^ _ 9_1_ fg T S\'- 
k 2ir g T \T c ) 



(dK)- 

bz e 



(1 



7k J 



7k) 

The latter equation is the first non vanishing correc- 
tion due to AL channel. Such a correction is always pos- 
itive, and it depends, as in the MT, on the lattice struc- 
ture factor 7k, but it does not vanishes in the regime 
T — T c ^> grS. This is a good feature of the system, 
since far from T c , the dynamical contribution plays an 
important role, and in this region, one has to compare 
it with DOS one, as discussed in the following section. 
Here, we just observe that since the corrections, Eqs. 
(|16|l . H21|) and (|30(l . have different signs, nonmonotonic 
behaviour in the total correction is expected, depending 
on the ratio gxS/T c . 



IV. DISCUSSION 

As we have seen, the total superconducting fluctuation 
correction to the thermal conductivity close to critical 
temperature is given by the following expression 



8 k 
kq 



TV 



4f / l./A'» 



:i-7k) 



3ir grS /i 
2 T c 



7k) - 1 



BZ 



z^(1-7k) 



(31) 

This correction has been obtained at all orders in the 
tunneling amplitude in the ladder approximation. Its be- 
haviour is plotted in Fig. 21 as a function of the reduced 
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FIG. 4: Total fluctuation correction to the thermal conductiv- 
ity for different values of grS/Tc for a two-dimensional system. 
A 1/e suppression is observed at high temperatures, with a 
sign depending on such ratio. At low temperatures, a finite 
correction, inversely proportional to the coordination num- 
ber 2, is reached at e = 0. In a finite interval of values of 
grS/Tc, a nonmonotonic behaviour of the correction is ob- 
served, where the correction is positive and increasing with 
decreasing temperature, reaches a maximum and then goes 
to a smaller (possibly negative) value at the critical temper- 
ature. 



temperature for the case of a two dimensional sample, 
and for different values of the ratio gr5/T c . We can rec- 
ognize two different regimes of temperatures: far from 
T c , e ^> gxS/Tc, and close to T c , e <C grS/T c . For a sake 
of simplicity, we will identify these two regimes as "high 
temperatures" and "low temperatures" , respectively. 

• High temperature regime e ^> gxS/T c . In this 
region, the electrons do not tunnel efficiently be- 
tween the grains and the system behaves almost 
as an ensemble of zero-dimensional systems. As 
a consequence, only the DOS and AL terms con- 
tribute significantly to the superconducting fluctu- 
ations; the correction to heat conductivity reads 



5k 3 6 1 

Kq IT 2 T c € 



3?r grS f 1 

2 T r 



1 



(32) 



This expression shows a 1/e singularity and it can 
have either positive or negative sign, depending on 
the ratio grS/Tc, we call 71 the value of the above- 
mentioned ratio solution of Eq. I|32fl . I n the ab- 
sence of renormalization due to tunnelling, the cor- 
rection is negative and corresponds to the typical 
singularity of the quasi-zero-dimensional density of 
state. On the other hand, increasing the barrier 
transparency grd, the correction grows due to the 
presence of the direct channel, i.e., the AL term, 
which becomes more and more important, until the 
correction itself vanishes at 71, after which it be- 
comes positive. A direct comparison with the be- 
haviour of the electrical conductivity 10 shows that, 
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already at this level, there is a positive violation of 
the Wiedemann- Franz law, being 



SL 



6k 5a 

3 9 g T Sz + l 7C(3) 



2tt T r 



• Low temperature regime e <§C gxS/T c . Here the 
tunneling is effective and there is a crossover to the 
typical behaviour of a homogeneous system, as T — > 
T c , from the point of view of the fluctuating Cooper 
pairs. Physically, the bulk behaviour is recovered, 
and one gets a non divergent (though nonanalytic) 
correction even at e = 0, where it equals 



5K(e = 0) _ 3 1 {3irg T 5 



Kq 



2 Tr 



- 1 



(34) 



The latter equation gives the saturation value in 
any dimension; it is also evident the 1/gr order of 
the perturbation theory. Again, the value of the 
constant can be either negative or positive. The 
correction vanishes at a value gxS/T c = 72 which is 
independent on the dimensionality and larger than 
7x- In the interval 71 < grS/Tc < 72, it has a non- 
monotonic behaviour, being positive and increasing 
for high temperatures and negative for low temper- 
atures. Such a behaviour has been represented, for 
the case of d = 2, in FigQJ The deviation from 
the Wiedemann-Franz law in the low temperature 
region is much more evident than in the high tem- 
perature one, because of the pronounced singular 
behaviour of the electrical conductivity close to the 
critical temperatureiS. 



V. CONCLUSIONS 

We have calculated the superconducting fluctuation 
corrections to heat conductivity. In the region of tem- 
peratures T — T c 3> 5t<5, a strong singular correction is 
found, reported in Eq. (|3*2*|) , corresponding to the sum of 
the DOS renormalization and the AL contribution in a 
quasi-zero-dimensional system. Moving closer to the crit- 
ical temperature, when T — T c -C gT&, the divergent be- 
haviour of the DOS term is cut off by the MT correction, 
which has opposite sign, while the AL term regularizes 
by itself to a finite value; this regularization signals the 
fact that the system undergoes a crossover to the homo- 
geneous limit. A nondivergent behaviour is found at the 
critical temperature, in agreement with previous calcu- 
lation in homogeneous superconductor o 6 ! 18 . The energy 
scale that separates the two regions, grS, can be recog- 
nized as the inverse tunneling time for a single electron 21 . 
As a final remark, we want to note that the ratio zgx5 /T c 
appears as the coefficient of the K— dependent term in 
the superconducting fluctuation propagator, Eq. <|11[) : 



from the standard theory of the superconducting fluc- 
tuations, the coefficient of K 2 in the propagator is ac- 
tually the superconducting coherence lengthi; we can 
therefore define an "effective tunneling superconducting 
coherence length" as £q T ^ = ayj grS/T c . From this defi- 

(T) 

nition, we can see that, if Q <§; a, the grains are strictly 
zero-dimensional at high-temperature and the correction 
to the thermal conductivity is always negative, while if 

(T) 

£ S> a, the direct channel of the superconducting cor- 
relations is strong enough to change sign to such correc- 
tion. 
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APPENDIX A: MICROSCOPIC DERIVATION OF 
FLUCTUATION PROPAGATOR 

Here we report a short description of the derivation 
of Ea. l|llll . evaluated in RefiiS, to remind the reader the 
main steps and the main assumptions of the calculation. 
We start from the expression of the partition function in 
the interaction representation 



Z = 



Tr exp ( - I H(r)dT 




Tr < exp 



xT T exp — 



H (T)dT 



H p {t) + H T 



dr 



(Al) 



We decouple the electronic fields in Hp by means of 
Hubbard-Stratonovich transformation, introducing the 
order parameter field A; because of our assumption, 
Et 3> A, the grains can be considered strictly zero di- 
mensional and we can neglect the spatial coordinate de- 
pendence in the field Aj in Eq. i|Al|l . We now expand 
over the field A,; the expansion is justified by our as- 
sumption to be close but above to the critical tempera- 
ture where the mean field (BCS) value of order parameter 
is still zero; moreover, we have to expand the action to 
the second order in t, too; this expansion is justified in 
the region 2 ^ t <C 1/t <C Et- We obtain two different 
contributions to the action: the first one is the typical 
action of superconducting fluctuations; the other one is 
the tunneling correction: ^off = S® s + Sl s . The first term 



V 



i - Attv f Tt ^2 ^ (e„, e M _ n ) 



2e ll >n„ 



(A2) 



8 



a 



b — c 

wvXKX 



Brillouin zone of reciprocal grain lattice. As it has been 
mentioned, the zero-dimensional character of the grain 
resides in the independence of the action on coordinates 
inside each grain. 

The tunncling-dcpendent part of the action is calcu- 
lated starting from diagrams in Fig0 they represent the 
first non-vanishing correction to fluctuation propagator 
due to tunneling. Their reexponentiation corresponds to 
the sum of the ladder series of tunneling and pairing in- 
teraction as reported in FigOO The calculation of diagram 
a, FigEfa), gives the contribution due to the possibility 
of tunneling of both electrons during the lifetime of the 
fluctuating Cooper pair, i.e. the Ginzburg-Landau time 
tql = 7r/8 (T — T c ); it is equal to 



FIG. 5: Top: The total tunneling correction to the fluctua- 
tion propagator self-energy is reported. The upper diagram is 
related to the possibility of tunneling of both electrons form- 
ing the fluctuating Cooper pair during the lifetime tgl of the 
Cooper pair itself. The other two diagrams consider the renor- 
malization of the intragrain fluctuation propagator. Bottom: 
The ladder series for the fluctuation propagator in the pres- 
ence of tunneling is reported. The crosses are BCS electron- 
electron interaction. 



s, 



t,(a) 



zgr V 7k |A K (fyiOT 



(A4) 



where, as mentioned, z is the number of nearest neigh- 
bors. 

Figures Efb) and^c), give an identical contribution, 
which is related to the probability that a single electron, 
participating in the fluctuating Cooper pair, undergoes a 
double tunneling, back and forth, during the Ginzburg- 
Landau time. Such a contribution reads 



0^, always appears as the combination of two fermionic 
Matsubara frequencies and it is therefore a bosonic one, 
as it should be. The sum over the fermionic frequen- 
cies in Eq. (| A2|) is logarithmically divergent and must 
be cut off at Debye's frequencji; using the definition of 
superconducting critical temperature, one obtains 

T 



in, 



In h lb - + J — ! =r 

T, V 2 AnT n 



where ip{x) is the digamma function, defined as the loga- 
rithmic derivative of gamma functioniS. Close to criti- 
cal temperature, T ~ T c , as already mentioned, the main 
contribution to singular behaviour comes from " classical" 
frequencies, |fi^| <C T c . Then, we can expand the ip func- 
tion in the small parameter \Q^\/T C : 



S eS =-*9T 



E I a k(O m )| 2 . 



(A5) 



The final result for fluctuation propagator at every order 
in tunneling in the ladder approximation is 



A K 



1 



V F In 



8T C 



zf^(l 



7k) 



(A6) 



Finally, we notice that the classical limit (f2 M = 0) for 
the fluctuation propagator Eq. (|A6|) can be obtained from 
a straightforward generalization of the Ginzburg-Landau 
functional for granular metals 

t m = E A dr ) A dr ') f°*< ( r ) *i ( r ') m ( r - r ') 



(A7) 



qO 



off 



- VF v S 



8T C 



|Ak (f2/i)| . (A3) 



In the last expression, for later convenience, we consid- 
ered the lattice Fourier transform: K belongs to the first 



where the parameter a is given by (l/4m£ 2 ) ln(T/T c ), 
where m is the electron mass, while the so-called Joseph- 
son parameter keep track of the tunneling effect: J — 
(\/Am^ 2 )(zg T 5/T c ). See also Ref. for the region of 
applicability of the theory reported above. 
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